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We derive a generalization of the Wiener-Khinchin theorem for nonstationary processes by intro¬ 
ducing a time-dependent spectral density that is related to the time-averaged power. We use the 
nonstationary theorem to investigate aging processes with asymptotically scale-invariant correla¬ 
tion functions. As an application, we analyze the power spectrum of three paradigmatic models of 
anomalous diffusion: scaled Brownian motion, fractional Brownian motion and diffusion in a loga¬ 
rithmic potential. We moreover elucidate how the nonstationarity of generic subdiffusive processes 
is related to the infrared catastrophe of l//-noise. 
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The Wiener-Khinchin theorem is a fundamental re¬ 
sult of the theory of stochastic processes. In its sim¬ 
plest form, it states that the autocorrelation function, 
C{t) = {x{t T)x{t)), of a stationary random signal x{t) 
is equal to the Fourier transform of the power spectral 
density, S{lu) = (|i(a;)p) /T [I112]- Here x{uj) denotes 
the Fourier transform of x{t), T the total measurement 
time and the ensemble average (...) is taken over many 
realizations of the signal. The Wiener-Khinchin theo¬ 
rem provides a simple relationship between time and fre¬ 
quency representations of a fluctuating process. Its great 
practical importance comes from the fact that the auto¬ 
correlation function may be directly determined from the 
measured power spectral density, and vice versa. Over 
the years, it has become an indispensable tool in sig¬ 
nal and communication theory [5], electric engineering 
[1], the theory of Brownian motion [3, classical [5] and 
quantum [7] optics, to name a few. 

The Wiener-Khinchin theorem is only applicable to 
weakly stationary processes that are characterized by 
a time-independent mean {x{t)) and a correlation func¬ 
tion C[t) that solely depends on the difference r of its 
time arguments [S]. It fails for nonstationary processes 
with an autocorrelation function, C(t, t) = {x{t-\-T)x{t)), 
that depends explicitly on time t and hence exhibits ag¬ 
ing. Further, the power spectral density, S{uj,t) —now 
an explicit function of both frequency and time—is not 
uniquely defined for nonstationary noisy signals [S]. A 
commonly used generalization, the Wigner-Ville function 
mm, given in Eq. Q below, is related to the instan¬ 
taneous power; it has the disadvantage of not being a 
true spectral density as it can take on negative values 
[5]. We here address these critical issues by first in¬ 
troducing a spectral density that is related to the time- 
averaged power. Contrary to the Wigner-Ville function, 
it is strictly positive and therefore a proper spectral den¬ 
sity. We employ this quantity to derive an extension 
of the Wiener-Khinchin theorem that is valid for finite¬ 
time, nonstationary signals. We apply this generalization 
to aging processes that are characterized by a correla¬ 
tion function of the form, C{T,t) ~ Cf^(j){T/1) , where 4> 
is an arbitrary scaling function. Correlation functions 


of this type describe scale-invariant dynamics that do 
not possess a distinctive time scale, contrary to expo¬ 
nentially correlated processes. They occur in a wide 
range of physical m, chemical [12] and biological [13] 
systems and have also found applications in finance m 
and the social sciences m- We illustrate our results by 
analyzing three important anomalous diffusion processes 
[16] : scaled Brownian motion, fractional Brownian mo¬ 
tion and diffusion in a logarithmic potential. We obtain 
for each of them analytical expressions for the nonsta¬ 
tionary spectral density in the low/high frequency do¬ 
main, and demonstrate the usefulness of the generalized 
Wiener-Khinchin theorem as a tool to perform spectral 
analysis in the aging regime. We finally discuss the gen¬ 
eral connection between subdiffusion and 1// noise [T7] . 

Nonstationary Wiener-Khinchin theorem. Let us be¬ 
gin by defining our notation. We consider a continu¬ 
ous random signal x(t) which is measured on the time 
interval t € As customary, we introduce the 

truncated Fourier transform of its time-shifted version 
vie) = xi9 + T/2) on 6» G [-T/2, T/2] [3]: 

rT/2 

yiu;,T)= dde-“^2/(0). (1) 

J-T/2 

Equation Q reduces to the usual Fourier transform in 
the limit T to infinity. The spectral density S{uj) is 
defined through the ensemble average of the modulus 
squared of y{uj,T) [3], 

^(w) = lim ;^(|y(a;,T)p). (2) 

T —>00 1 

Note that Eq. (H has to be interpreted in the distri¬ 
butional sense for nonintegrable functions. In the fol¬ 
lowing, we assume that the stochastic process x{t) is 
real. As a result, the power spectrum S{uj) and the 
stationary correlation function C'(t) are even functions. 
Physically, the spectral density measures the power con¬ 
tained in the signal at a certain frequency [3]. More 
precisely, duj S(lo)/t: corresponds to the power con¬ 
tained in the frequency interval [wi, UJ 2 ] and the total 
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power {x^) is obtained by integrating over all frequencies, 
{x^) = f^duj S(uj)/(27r). According to the Wiener- 
Khinchin theorem S(u!) and C(t) are related via [HE], 

/ oo ^OO 

dre““'^C'(T) = 2 / dr cos(wt)C(t). (3) 

-OO J 0 

Equations (§ and (§ may be regarded as two definitions 
of the power spectrum. They are fully equivalent for 
stationary processes |1|. 

For nonstationary random signals, Eqs. ([^ and ([^ are 
no longer equivalent and both expressions may be used to 
define a nonstationary generalization of the power spec¬ 
trum. Starting from Eq. ([^, the Wigner-Ville spectrum 
S'wv(w,t) is constructed as the Fourier transform of the 
symmetrized correlation function [18j . 

S'wv(w,t)=y At e~^'^'^{x{t+^x{t-^Y (4) 

Equation @ is simply related to the nonstationary corre¬ 
lation function, C{T,t) = /^(dw/27r)e*‘^'^5'wv(<j-', <), thus 
preserving the usual form of the Wiener-Khinchin theo¬ 
rem. However, it is not a proper spectral density since it 
can take on negative values. Another drawback is that 
it depends on the infinitely extended past and future of 
a process. This does not pose a problem if the investi¬ 
gated signal is non-zero only in a finite time interval. But 
the truncation of the correlation function of an extended 
signal which results from a finite measurement time will 
lead to deviations from the Wigner-Ville spectrum that 
are not taken into account in Eq. Q. 

We here follow a different path and interpret Eq. ([^ 
as the definition of the spectral density. Thus, for a finite 
measurement time T, we introduce the quantity, 

5ta(a;,r) = i(|y(a;,T)|2) 

/.T/2 ^T/2 

= -/ d0' d0" 

J J-T/2 J-T/2 

(5) 


This definition ensures that 5'ta(w, t) is positive and hence 
a true spectral density, contrary to the Wigner-Ville spec¬ 
trum Q. For stationary systems, we have, 

S'(a;) = lim S'ta(w, T) = 5'wv(w, 0), (6) 


so that all three definitions of the spectral density coin¬ 
cide in the long-time limit. For nonstationary processes, 
5 ( 0 ;) does not exist, whereas 5'ta(w,T) and S'wv(w,t) are 
generally distinct. The Wigner-Ville spectrum Syjv{uj,t) 
is related to the instantaneous power, while 5'ta(w,T) cor¬ 
responds to the time-averaged (ta) power. 


1 

27r 

1 

27T 



doj S'wv(w,i) = {x^{t)), 

dw Sta{uj,T) = ^ [ dt' {x^{t')). 

^ JO 


(7) 

( 8 ) 


We note that Eqs. Q and (|^ differ for nonstation¬ 
ary signals. Equation (|^ provides a clear physical 
interpretation of 5'ta(w,T): f^^duj 5'ta(w,r)/7r is the 
time-averaged power contained in the frequency interval 
[ 0 : 1 , 102 ], in complete analogy to the stationary case. 

From the definition ([^, we find that the time-averaged 
spectral density S'ta(w,T) can be expressed in terms of 
the correlation function C{T,t) as. 


SU^,T) = ^ f dt" [ dt'cos{iot')C{t\t"-t'), (9) 

where we have used that C{t, t) = C{—t, t + t) so that 
both arguments of the correlation function are positive. 
On the other hand, the autocorrelation function is related 
the spectral density S'ta(w,t) via. 


C{t,T-t) 


1 + T 


d_- 

&f. 



dw cos(a;r)S'ta(w, T), 

( 10 ) 


for T > T. Equations 0 and ( |I0| ) constitute a 
nonstationary, finite-time generalization of the Wiener- 
Khinchin theorem ([^. It is free from the problems faced 
by the Wigner-Ville function Q ; the price to pay is that 
the familiar symmetry between time-frequency represen¬ 
tations known from the Fourier transform is lost. In the 
limit of stationary signals, Eqs. 0 and ( [To| ) reduce to 
the usual Wiener-Khinchin theorem 0. 

Application to scale-invariant processes. In order to 
illustrate the usefulness of the nonstationary Wiener- 
Khinchin relation ([^-(10), we next consider aging pro¬ 
cesses with a correlation function of the scaling form m, 




( 11 ) 


Here C > 0 is a constant, a > —1 the scaling exponent 
and 4>{u) an arbitrary scaling function. Rescaling of the 


time variable only changes the prefactor of Eq. (11) and 
not its qualitative behavior. This kind of autocorrelation 
function is understood to describe the long-time behavior 
of a system, where both the age t and the time lag r 
are large compared to the system’s intrinsic time scales. 
In the absence of a characteristic correlation time, the 
only relevant time scale is the age t of the system, which 
also governs the decay of the correlations as a function 
of the time lag r. Scaling correlation functions of this 
type appear in a large number of anomalous diffusion 
problems as we will discuss in the next section. 

We may now use Eq. 0 to compute the time-averaged 
spectral density 5'ta(w, T) for the scaling correlation func¬ 
tion 0- We immediately find, 

5ta(w,r)~2CT“+ix(wr) with 

pi nOO 

x{z) = J du cos(u 2 ;) J du (12) 
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Remarkably, the time-averaged spectral density has a 
scaling form, similar to the autocorrelation function 0, 
with a scaling function y;(a;T). Using Eq. (101, we further 
obtain the corresponding inverse relation, 

1 

(j){u) = - dz cosizu)[{a+ 2)x{z) + zx'iz)]. (13) 
Jo 


The nonstationary Wiener-Khinchin theorem (§-@ 
thus leads to a direct relationship between the time and 
frequency scaling functions (j){u) and x(z). While the be¬ 
havior of x{z) depends in a nontrivial manner on 4'{u) 
and the exponent a, the mere existence of the scaling 


form (12) already implies a number of generic proper¬ 
ties for the spectral density We first note 

that the scaling function xi^T) depends on the prod¬ 
uct of frequency and measurement time. This means 
that low and high frequency limits of Sta{u},T) are in¬ 
timately connected to the measurement time T. For 
low frequencies such that uiT 1 (or, equivalently. 


short measurement times for a fixed frequency), Eq. (12) 


may be easily evaluated since x(0) is just a constant, 
and we obtain S'ta(w,r) ~ 2CT“+^x(0)- In the low- 
frequency limit, the spectral density hence increases al¬ 
gebraically with the measurement time and is indepen¬ 
dent of frequency to leading order. This reflects the fact 
that frequency components, whose period is much longer 
than the measurement time, will be essentially constant 
over the measurement interval. For these frequencies, 
the truncated Fourier transform Eq. 0 reduces to the 
time average x{T) = f^dt x{t)/T of the process x(t). 
Since by definition, S'ta(0,T) = T{x^{T)), we hnd that 
{x^{T)) ~ 2Cx(0)T““^. This implies that processes de¬ 
scribed by the scaling correlation are markedly dif¬ 
ferent depending on the value of the a. For a < 0, pro¬ 
cesses are ergodic, as the variance of the time average 
tends to zero, whereas for a > 0 it increases with time 
and processes are generally nonergodic [20] . 

In the high-frequency regime, ujT 3> 1 (or large mea¬ 
surement times for a given frequency), we may explicitly 
compute the asymptotic expansion of the function x(z) 
for z > 1: 


X{z) ax 


expansion of the spectral scaling function will be propor¬ 
tional to as the first term in Eq. (14) vanishes for 


A = 0. As a result, the nonstationary spectral density, 
Sta,{uj,T) oc exhibits a uj~'^ frequency depen¬ 

dence like a normal diffusive process, which is recovered 
for a = 1. By contrast, if the expansion of (/'(u) con¬ 
tains a nonanalytic term with A ^ No, the corresponding 
term in the expansion of xi^) proportional to z~^~^ will 
be non-zero. In this case, to leading order, the spectral 
density behaves as S'ta(w,r) oc This kind 

of frequency dependence is generally referred to as 1//- 
noise m- We will elaborate on this connection in more 
detail below. 

Examples. Let us apply the above results to three 
paradigmatic models of anomalous diffusion. These pro¬ 
cesses are characterized by a nonlinear increase of the 
mean-square displacement, (x^(t)) = 2Dat°‘, where 
is a generalized diffusion coefficient. Subdiffusion is ob¬ 
tained for 0 < a < 1 and superdiffusion for a > 1 [TH] . 

(a) Scaled Brownian motion (SBM). SBM is maybe the 
simplest way to construct a stochastic process that ex¬ 
hibits anomalous diffusion. It is defined as x^^^{t) = 
a;(t“), where x{t) is a Brownian motion [5^ 123] . Its au¬ 
tocorrelation function reads for r > 0, 

C7®™(t, t) = 2Li„min(f“, {t + t)“) = 2D^r, (15) 


It is of the scaling form ( |ll| ) with C = 20^ and 
^SBM(.^) _ ^ SBM, like Brownian motion, is Marko¬ 
vian, however, its increments are nonstationary. 

(b) Fractional Brownian motion (FBM). FBM is a non- 
Markovian generalization of Brownian motion. One usu¬ 
ally distinguishes between Riemann-Liouville (RL) FBM 
with nonstationary increments [24] . 


„RL-FBM 


{t)= f dx(t'){t-t')^-\ (16) 

Jo 


and Mandelbrot-van Ness (MN) FBM with stationary 
increments 1251. 


„MN-FBM 


{t) = f dx{t')(t — t') ^ ^ 
Jo 


(17) 



+r 

a-A-bl V 2 / 

J — oo 




-l- r(A -l- 2) cos 




_l_ 0{z~^~^) (14) corresponding scaling autocorrelation functions are 


given by Eq. (11) with 


provided that the scaling function asymptotically be¬ 
haves as (j){u) ~ axu^ with A > — 1 for u ^ 1 (see details 
in the Supplementary Material m)- We observe that the 
frequency dependence of the spectral density S'ta(w, T) in 
the long-time limit depends on the small argument ex¬ 
pansion of the scaling function (/)(u) and, hence, on the 
long-time behavior of the correlation function for t 3> t. 
If the scaling function is analytic for small arguments, i.e. 
(j){u) ~ oq -f aiu + 0{u^), then the leading term in the 


/■l C.-1 

^RL-FBM(y^ ^ ^ / dww—{w + u) = (18) 

^MN-FBM(y) ^ (19) 

SBM and FBM motion play an important role in describ¬ 
ing anomalous diffusion in biological living cells pSl - BU] . 

(c) Dijfusion in a logarithmic potential (LOG). A Brow¬ 
nian particle diffusing in an asymptotically logarithmic 
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potential of the form U{x) ~ ?7oln|a:/a|, for x ^ a, 
with free diffusion coefficient D* exhibits subdiffusive be¬ 
havior for 1 < Uo/D* < 3, with a diffusion exponent 
a = 3/2 — Uq/2D* [30]. The autocorrelation function is 
also of the scaling form (11) with a scaling function m, 




/Tra 


r(3-a) 


( 20 ) 


^OO o V 

J dw w^e~'^ iFiy-,3 — a,w^jT(2 — a,w^u). 


Diffusion in logarithmic potentials appears in a great va¬ 
riety of problems [32] . ranging from DNA denaturization 
[34], systems with long-range interactions [33] and diffu¬ 
sion of atoms in optical lattices [35] . 


In the following, we compute the spectral density (§ 
for the three processes by applying the nonstationary 
Wiener-Khinchin theorem (12). For SBM, the scaling 
function is unity and thus analytic. By Eq. (14), we then 
have x(z) cx for large z (see Fig. [^. Consequently, 
the spectral density is for ujT 1, 


cSBM 

‘^ta 


{oj,T) ~ 4DaW-2T“-\ 


( 21 ) 



FIG. 1. (color online) Frequency scaling function x(z) (121 
for four models of anomalous diffusion: scaled Brownian 
motion (SBM), Riemann-Liouville (RL) and Mandelbrot-van 
Ness (MN) fractional Brownian motion (FBM) and diffusion 
in a logarithmic potential (LOG) for the scaling exponent 
a = 0.3. For identical mean-square displacements in all cases, 
the asymptotic dependence of x(^) for large 2 (dashed lines) 
is different, owing to the nonanalytic behavior of the time 
scaling function (jiiy) for both types of FBM and LOG. The 
results shown here are exact results obtained by numerical 
integration of Eq. (121, the asymptotic expansion for large 2 
yields Eqs. (211 through (24l. 


with an algebraic dependence on the measurement time 
T. If the process is subdiffusive, 0 < a < 1, the magni¬ 
tude of the high-frequency spectrum decreases with time, 
as shown in Fig. [^ The rescaling t ^ in SBM can thus 
be understood as a continuous slowing down of the pro¬ 
cess, which reduces the prevalence of high-frequency com¬ 
ponents in the spectrum. In the superdiffusive regime, 
1 < a < 2, the process speeds up and the magnitude 
of the high-frequency spectrum increases with time. In 
both cases, the dependence on frequency is precisely the 
same as for usual Brownian motion. 

On the other hand, for both kinds of FBM and for 
the logarithmic potential, we have (/(u) ~ oq -b aaU°‘ + 
aiu (the explicit expressions for the coefficients are given 
in the Supplementary Material m), and hence x{^) oc 
2 “““i for a < 1, see Fig. The time-averaged spectral 
density in the subdiffusive regime, 0 < a < 1, and for 
long measurement times, ojT 3> 1, is to leading order 
given by. 


allows to evaluate finite-time corrections to leading order, 
c-FBM.LOG/ rr\ nr -a-1 

(w,T)~2Cw aar(Q:-b l)sm J 
-b aar(Q! -b 1) cos 

+ (ao-^)(wT)“-i , (23) 


These corrections have two different origins: The domi¬ 
nant terms proportional to oq and oi stem from the non- 
stationarity of the correlation function. The correction 
term proportional to Oq, on the contrary, appears even 
if the system is perfectly stationary (e.g. for a < 0) and 
is due to the finite measurement time. In the superdiffu¬ 
sive regime, 1 < a < 2, the spectral density of FBM also 
displays an algebraic time-dependence to leading order, 

S'f™(a;,r) ~ 2c(ao - (24) 


r.FBM,LOG/ N nr v/ i -a-1 

Ata (w) ~ 2Caar(a-b l)sin —Jw 


( 22 ) 


Interestingly, even though the process is nonstationary, 
the spectral density is asymptotically independent of the 
measurement time, see Fig. [^ This shows that a sta¬ 
tionary spectral density does not necessarily imply a sta¬ 
tionary process, the nonstationarity being in this case 


encoded in the low-frequency cutoff. Equation (14) also 


The overall magnitude of the spectral density increases 
with the measurement time T. This is similar to the 


result (21) for superdiffusive SBM. 


Relation to Ijf-noise. As we have seen in Eq. 
both kinds of FBM and the diffusion in a logarithmic po¬ 
tential lead to a spectral density of the form 5'ta(w) oc 
. Generally, a spectrum with S{oj) oc u}~^ with 
0 < /3 < 2 is termed 1//-noise [T7]. This kind of spec¬ 
trum was first observed in flicker noise in triodes |36| , and 
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FIG. 2. (color online) Time-averaged spectral density 
Sta.{iJ,T) (121 for scaled Brownian motion (SBM) and 
Mandelbrot-van Ness fractional Brownian motion (MN-FBM) 
for a = 0.3 and different measurement times T. The spectral 
density converges to a stationary expression in the long-time 
limit for FBM, whereas for SBM, the magnitude of the high- 
frequency spectrum decreases with time. 


has since been found in a wide range of other systems 
in physics m, biology [3S], geology [SU] and phonology 
m- The scaling Wiener-Khinchin relation (U^ shows 
that l//-noise naturally occurs for processes with scale- 
invariant correlations. The main requirement is that the 
correlation function is nonanalytic in the limit t ^ t. 
Out of the three examples we investigated above, this is 
the case for FBM and the diffusion in a logarithmic po¬ 
tential, with the former being non-Markovian but linear, 
and the latter Markovian but nonlinear. For a < 0, these 
processes are stationary and the power spectrum, which 
is proportional to is integrable at low frequen¬ 

cies, yielding a finite stationary power. For 0 < a < 1, 
by contrast, the stationary spectrum is nonintegrable at 
low frequencies (infrared catastrophe), formally leading 
to infinite power. This has caused some discussion on the 
physical nature of this divergence m- Using the scaling 
Wiener-Khinchin relation (12), this apparent contradic¬ 
tion is easily resolved, since for a diffusive process, the 
power actually diverges in the infinite time limit. For 
any hnite time, there is a low-frequency cutoff on the 
spectrum that depends on the measurement time, which 
ensures finite power. As the measurement time increases, 
this cutoff moves to lower and lower frequencies, and 
the process approaches 1//-noise at all frequencies. We 
note that a similar behavior was recently observed for an¬ 
other scale invariant (superdiffusive) stochastic process, 
the Levy walk, in application to blinking quantum dots 
dH 133] . Our approach allows to generalize these results 
to generic subdiffusive random processes. 

Conclusion. The generalized Wiener-Khinchin theo¬ 
rem ([^-(10) extends the notion of a power spectrum to 


nonstationary stochastic processes in a natural way and 
establishes a direct connection to their autocorrelation 
functions. By studying three ubiquitous scale-invariant 
models of anomalous diffusion, we have demonstrated its 
usefulness in analyzing the properties of aging processes. 
Our results thus provide an extension of the range of 
applicability of spectral analysis—in terms of a proper 
(positive) spectral density—to the nonstationary regime. 
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SUPPLEMENTAL MATERIAL 


Large argument expansion of x(^) 


In order to obtain the long-time limit of the spectral density, we need to know the asymptotic behavior of the 
function x{z) for large arguments. From Eq. (12) of the main text, x(z) is defined as, 


p1 poo 

x{z) = J dv cos{vz) J drt 


(S25) 


We define the auxiliary function. 


/ OO 

du m“““^(/)(m). 


(S26) 


For the scaling function 4>{u), we assume the asymptotic behavior, 


(/)(m) ~ 


a\u^ for M <C 1 
for 1, 


(S27) 


with A > — 1 and < a + l- Using this, it is easy to show that. 


f{v) ~ 


a\ 




a+l-A 


for V <C 1 




— _ H _ (1 — 


for 1 — u <C 1. 


(S28) 


We now split the integral into three parts. 


re rl-e /•! 1 


10 Je Jl-e 


du cos{vz)f{v) 


(S29) 
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where e is small enough that the expansion (S28) is justified. Using this expansion, we find for large z, 
du cos(uz)/(u) ~ ua —^^ — r(A + 1) sin + 

^A+l , 


e'' sin(ze) ^ Ae'^ ^ cos(ze) 
z 


- r(A + 2)c» tu 

V 2 / z z^ 


du cos(uz)/(u) ~ 6 


/l-£ 


1 


— r(Q; + 2 — cos ^z 


r(a — /r) 


^ a + 1 — fi 

(a + 1 — cos(z(l — e)) 






ga+l-Msin(z(l - e)) 


+ 




a + 1 — /r 


r(a + 3 — /r) sin ^2 


r(a - 


u—a — 3 


e“+2-^sin(z(l - e)) 


+ 


(a + 2 — ^)e“ cos(z(l — e)) 


(S30) 


In the center part of the integral, we may integrate by parts. 


[ du cos(uz)/(u) = 

J e ^ 


cos(uz) , 


/ 


■C>(z 


-3\ 


(S31) 


Again using Eq. (S28) for the boundary terms, we see that all terms depending on e cancel and we are left with, 
X(z) ~ ax 


„ +1 _ A ™ 


+ 


r(A + 1) . -x -1 

b, ' 

a + 1 — /i 


— r(Q! + 2 — ^) cos ^2 


7r(Q; — /i) 




r(Qf — /i) 


fj,—oc—3 


. (S32) 


For ^ < a — A, which is the case for all the examples discussed in the main text, the leading order expression is given 
by the part proportional to oa, which is precisely Eq. (14) of the main text. 


Asymptotic expressions for (f>{u) for the example systems 



{>4 




SBM 

0 

1 

0 

1 

RL-FBM 

0, Q, 1 , ... 

1 . V^r ((a + l)/2) _ CK. 

Q —1. a—3. ck—5 . 

2a , Q:(q: —1). 0 ( 0 : —1)(ck —3) . 

■*-’ r{a/2) sin(7rci/2) ’ 2 ’ * ’ ' 

2 , 2 ’ 2 ’ ■ • • 

a + 15 ck+3 ’ 4(q;+5) ’ 

MN-FBM 

0, Qi, 1, ... 

1* —1- -• 

j., j., ^ . 

0; 0 — 1; a — 2; ... 

1 , a . Q:(q: —1) , 

2 ’ 2 ’ 4 ! • • • 

LOG 

0; a; ... 

x/t: _ 7rr(2 —Q:)r(l —q:) ^ 

UL 2 ’ 2 ’ ■ ■ ■ 

V^ar(7/2-Q:) . x/^Q:(3-2a)r(9/2-a) . 

2r(3-ct)’ 4r2(3/2-a) ’ ••• 

3r(3-a) ’ 10(a-3)r(3-a) ’ ' ' ' 


In the main text, we discuss the spectral density for scaled Brownian motion (SBM), Riemann-Liouville (RL- 
FBM) and Mandelbrot-van Ness (MN-FBM) fractional Brownian motion and diffusion in a logarithmic potential 
(LOG). The table lists the explicit expressions for the coefficients and exponents of the small argument {(j){u) ~ 
OAiU^i + ax^u^'^ + ...) and large argument ((^(u) — + ...) expansions of the corresponding scaling 

functions. The first column lists the stochastic model, the second column the first few exponents Xi in the small- 
argument expansion in ascending order (for 0 < a < 1). The third column lists the coefficients corresponding to these 
exponents. The fourth and fifth column are the same, but for the large-argument expansion. So, for example, the 
first three terms in the small argument expansion for the scaling function for MN-FBM motion are, 

(/)(u) ~ 1 — u“ -I- -I- 0{u^). (S33) 

Note that each of the exponents Xi and /ii individually yields a contribution of the type (S8) to the expansion of the 
frequency scaling function x(z). 






























































